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Abstract 



This work present a new class of variational wave functions for fermi systems 
in any dimension. These wave functions introduce correlations between Cooper 
pairs in different momentum states and the relevant correlations can be computed 
analytically. At half filling we have a ground state with critical superconducting 
correlations, that causes negligible increase of the kinetic energy. We find large 
enhancements in a Cooper pair correlation function caused purely by the interplay 
between the uncertainty principle, repulsion and the proximity of half filling. This 
is surprising since there is no accompanying signature in usual charge and spin 
response functions, and typifies a novel kind of many body cooperative behaviour. 



Introduction: Variational wave functions have played an important part in our un- 
derstanding of several phenomena in quantum manybody systems. The role of the BCS 

wave function in superconductivity, Jastrow-type functions for superfluidity 0, the 
Laughlin function in the Fractional Quantum Hall problem ||], and the Gutzwiller wave- 
function in strongly correlated fermi systems have been recognized to be of fundamental 
importance. 

This work presents a new variational wave function and some variants of it, that 
promise to be of interest in the topical problem of superconductivity arising from repulsion. 
We have recently argued on the basis of certain inequalities ||, that projecting out s-wave 
Cooper pairs in a fermi system on a lattice would lead to enhanced extended s-wave type 
pairing fluctuations near half filling. These are expected to lead to superconductivity at 
precisely half filling for a class of Hubbard-type models. In brief the argument involves 
the recognition that s-wave and extended s-wave pairing are canonically conjugate in 
the sense of the uncertainty principle, which ultimately drives the instability towards 
superconductivity near half filling. This conjugacy arises since the s wave and the extended 
s wave Cooper order parameters B, A ( defined below) satisfy the commutation relation 
[B, = 2T, where the kinetic energy operator T plays a benign role, similar to that of the 
number operator N as in other familiar contexts, such as the conductivity sum rule on the 
lattice. From the uncertainty principle, it follows^] that the fluctuation < A^A + AA* > 
is bounded from below by 4 < T > 2 / < B^B + BB^ >. The special role played by 
half filling is due to the fact that the suppression of s-wave fluctuations is possible to a 
very high degree near half filling, the commutation relation [B, B'] = L — N in a sector 
< N >~ C , permits both < B^B > and < BB^ > to be simultaneously small. This 
"squeezing" results in large enhancements in the conjugate variables, remarkably enough, 
without an appreciable energy cost. This effect has been termed order by projection. 

While the arguments in || provide a novel direction, the absence of explicit solutions 
or of good wave functions has proved to be a hindrance in arriving at a thorough under- 
standing of the physics of these models. We provide in this work a wave function which 
is tractable enough so that calculations of all relevant expectation values can be done 
analytically, which by itself is rare in many body systems. Further the wave function 
catches the essence of the enhancements mentioned above, and at half filling is argued to 
be exact. While we know of no materials to which these models apply, we make bold to 
present these models and wave functions, since they address and provide what we believe 
to be amongst the first positive results |§ for the important theoretical question of the 
possibility of superconductivity without explicit attractive interactions- in fact arising out 
of pure repulsion. 

The model and the variational wavefunction: Let us write the Hamiltonian of 
the model considered in d dimensions as 

H = T + U S B ] B . (1) 

Here B = °k is the s-wave Cooper pair operator, and bk = c_fcjCfcf the pair operators, 
the kinetic energy T = J2 ^k^ka, with the band dispersion = —2 X^=i cos(k a ), although 
later we will consider a slightly more general form of the band dispersion in 2-dimensions 
with orthorhombic distortion. The U s term discourages on site s-wave correlations, and in 
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fact projects out s-wave order. In order to get a feeling for its effect, consider the estimate 
of the energy in a BCS state \^bcs >— Y\{uk + Vkb\)\vac >. Provided J2 u kVk = 0, we find 
Ebcs = S[2(efc — n)v\ + U s Vff\, subject to N = 2 v\. The role of U s is clearly to flatten 
out v k from its step function behaviour in the normal state, and hence a superconducting 
state with true Long Ranged Order (LRO) and an extensive energy shift arises, at least 
as a local variational minimum. In Ref||, a more general model is introduced including 
the above term as well as the Hubbard U term. The general results of f| imply for these 
models that away from half filling, the ground state energy density is unshifted by the 
U s term, while the extended s-wave correlations get a large enhancement at the expense 
of suppressing s-wave correlations. By continuity in density, we expect superconductivity 
at half filling. In this work we ignore the Hubbard U term for tractability, and explore 
the other aspects within the model of Eq ([[]), in particular the origin of the enhancements. 

The extended s-wave operator A = —2 J2 satisfies the commutation relations 
[-B,T] = —A and [B, A*] = 2T. The main variational wave function proposed here is 
written in terms of the free fermi wave function |$ >= ]TT | fc| <fcy b\\vac > as 

\* > = exp(-^ S) |$ >, (2) 

where the number of lattice sites is C. The prefactor generates several Cooper particle 
pairs and Cooper hole pairs in the fermi gas. To motivate this wavefunction, note that 
the anti Hermitean operator S = £)( e fci — ^Wk-fikz can be viewed as the commutator 
\ [T, B^B] , and hence the wave function may be viewed as a "rotation" about a direction 
orthogonal to the kinetic and potential energies. Such a strategy is familiar e.g. from 
diagonalizing a sum over two Pauli matrices, as well as quadratic forms in bosons. A 
similar rotation also gives the exact answer in a further simplified model ||. 

Calculation of expectation values: We now turn to the calculation of expectation 
values < Q >$=< ^e\Q\^e > of various operators Q in the above state. We begin by 
evaluating the derivative of the kinetic energy < T >g, 

-| < T > e = h< A^A > e -2 < &B + StC > e ] = 4a(0) - 2v{6), (3) 
do L. 

where C = J2 e t^k, thus defining the functions a, v. We immediately note a feature of 
considerable importance, namely that unless one of the correlation functions in the right 
hand side of Eq(||) possesses explicit Long Ranged Order (LRO), the change in the ( 
extensive) kinetic energy is not extensive, and may be neglected in the thermodynamic 
limit. We will make repeated use of this insight later. We also note that this includes the 
case where one of the correlators possesses Quasi LRO, i.e ~ with C < 1, indeed we 

find later that C = | does occur for the case of greatest interest, namely the half filled limit. 
We also need the expectation of operators T n = e k 4>k, with (fi k = n k ^ + n_ k ^ — 1, so that 
T\ = T. The equations are most compactly written in terms of the operators I n = cfflk 
such that Iq = B, Ii — —\A, I2 = C, and their correlators r n>m (9) = 7 < I^I m + I^In >e 
as 

-!<T n >*=2[T nil (0)-r n ,o(0)]. (4) 
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We note that in view of the above definition, LRO corresponds to T ~ 0(C), whereas 
in the absence of true LRO, T ~ o(C), and hence here the correlators < T n >g are 
insensitive to the rotation, at least to leading order in C. 
The correlators r„ jm satisfy exact equations 

4^(0) = -L < [T n+1 lU m + I%T m+1 - l\hT m - Tj\l m + h.c] > , (5) 

with initial conditions obtained by taking the expectation value in the free fermi ground 
state: they read r n m (0) = 2[e k +m f k ] k , with f k the usual ( non-interacting) fermi function 
and [g k ] k = \Y.k9k- 

The equations are handled assuming a fundamental factorization that arises in the 
thermodynamic limit. This factorization of hermitean global operators ( i.e. sums over 
all sites of local operators) and may be expressed as the statement that for any two such 
Qj with nonzero averages, we have < Q1Q2 >e /[< Qi >e< Q2 >$] ~ This is 

also equivalent to the statement that the connected part of the correlator is subleading 
in powers of C. The case needed by us corresponds to Q\ ~ £ T n m and Q 2 ~ T n . If this 
is assumed then we immediately see that Equations® factor out into linear equations 

with coefficients \i n = [e^(2fk — l)]k obtainable numerically, together with the initial 
conditions i\ m (0) = 2[e n k +m f k } k . 

We have checked this assumption to low orders by evaluating the correlation functions 
r nm for the most relevant cases (n,m) = (0,0), (1,1) out to 0{6 b ). The calculation 
proceeds by expanding the expectation value as a nested commutator with S, which is 
evaluated by Wick's theorem for ground state correlations, leading to n\ contractions that 
are further multiplied by powers of e. This results in sums over k of algebraic functions 
of the band energies and fermi functions. The series agrees with the one generated from 
the above factorization procedure exactly to that order, and provides non trivial support 
for it. 

We now consider the correlation functions of interest, and begin by denoting (3(9) = 
|r 0) o, Oi(6) = Ir^i, 7(6*) = r^o and finally v(9) = r 2 ,o- The equations of immediate 
concern are then 

2^(3 + (1 - p) 7 

2(1 - p)a + 2/i 2 /3. (7) 

We also note for later use that 4sV = [i\v + 2/z 3 /3 + (1 — p)i\ 2 — /U 2 7- We denote the 
particle density as p such that p < 1. The initial conditions are /3(0) = \p, a(0) = u(0) = 
d+ |/U 2 , 7(0) = p\. We now examine the behaviour of these coupled equations. Let us 
note that po = (p ~ 1) ; further p\ is negative and non vanishing near half filling, whereas 
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fj,2 also negative, is very small ~ 0((1 — p) 3 ) near half filling. In fact for all even n we 
must have \i n — > as we approach half filling, due to particle hole symmetry. If we neglect 
the cross couplings, then (3 falls off as 9 increases from while a increases. 

A neat way to solve these equations is suggested by inspection: we observe that 
in effect the Eqns(^) are generated by a linear relation 4sl n = p n +ih — Huh and the 
hermitean conjugates of these. The 9 dependence of all the operators are related thus 
to the fundamental ones Io,I±, and thus we find linear combinations that have simple 



evolution: Q± = —poh + (pi ± Aq)/q where A = yp? ~~ MoA*2- These obey uncoupled 
equations ^Q± = ±A Q±. 

We readily invert and find with £± = (1 ± I = ^{Q+ — Q-} and I\ = 
2(i-p) i£-Q+ + £+Q-}- The four basic correlators F ait(T2 (9) = < Q^Q^ >e can thus 
be found simply 'F ai ^ 2 (9) = F CTl)O . 2 (0) exp(Ao#(o"i + 02)). The initial conditions are given 
as 

F ai!<72 (0) = (l-p) 2 a (0) + (l-p){ /Ul + ^( ( 7 1 + a 2 )}7(0) + (p 1 + ( T 1 Ao)(/ii + a 2 AoMO). (8) 
We thus find the expectation values: 

<x{9) = 4(1 \ )2 [e^i,i(0) exp(2A g) + e + F-i,-M exp(-2A fl) + 2^-^1,-1(0)] 

(3(9) = -^[F hl (0) exp(2X 9) + F_ 1( _i(0) exp(-2A o 0) - 2F 1) _ 1 (0)] 

7(0) = 2Xq{ 1 _ p) [^1,1(0) exp(2A g) - £ + F_ w (0) exp(-2A g) + 2(^)F W (0)I9) 

These together with Eq(JJ) provide the solution to the variational problem since the vari- 
ational energy density is \i\ + U s /3(9), as well as that of computing the correlations. 

We show in Fig(l) the behaviour of the three functions a, (3, 7 for the case of 2-d at a 
typical value of the density p = .75. Note that (3 goes through a very shallow minimum 
with a value (3*(p) = (3(9*) ~ 0, and 7 also seems to go through zero at nearly the same 
value of the minimizing 9 = 9*. From the Equations^, (3* would be actually zero if 7 
vanishes at exactly 9*. To investigate this further, we plot in Fig(2) the minimum value 
(3*(p) as a function of density in 2-d, and the inset shows a similar plot for 1-d. We see that 
for general densities, although (3* is impressively small, it is non zero. From the general 
arguments of || we know that the exact ground state energy density is unshifted from 
the non interacting value, corresponding to (3 exact = 0, and so our wave function gives 
a very good approximation to the true energy for most densities, provided we restrict 
to U s ~ \pi\. In Fig(3) we show a*(p)/a(0), this is enhanced greatly, by an order of 
magnitude as far as 15% away from half filling. The lower bound = |( rfz^ + ^2) 
obtained in Ref|5|], also has a similar behaviour close to half filling. 

Solution at Half filling: The discussion of the solution above indicates that for 
half filling, i.e. p = 1, the wave function discussed here could be exact, since (3* goes to 
zero here. Our analysis of the orders of magnitude of the correlators was based so far 
on the implicit assumption that 9 ~ 0(1). From the solution it can be shown that for 
p < 1, the minimizing 9* ~ — log(l — p), and hence at p = 1, we must work with greater 
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precision to see if the assumptions made earlier are consistent. Firstly we assume that 
Eqns(^|) are still valid at half filling, requiring that lim^^ < T n > e /C is unchanged 
from its non-interacting value We verify at the end that this is still true and hence 
obtain a self consistent argument. With the assumption we set fi2n and hence the 
equations decouple and are integrated immediately giving a{9) = a(0) exp(— 2/ii6>) and 
(3{9) = (3(0) exp(2/i 1 6 l ). We also integrate the equation for v and find v[ff) = g 1 exp(/i!#) + 
gi exp(2/Z!#). The constants g\^g% involve the initial conditions, but are not important 
since we see that for large (positive) 9 we can safely neglect u(9) ( recall that \L\ < 0). 
We substitute for a into Eqn([|) and integrate over 9, to find the energy per site0 

e(6) = e non + ^ exp(2//!0) + exp(-2/ii#) + (negligible terms) (10) 

Z \fii\jL 

where e non = Hi is the ground state energy density of the free fermi gas, and the negligible 
terms are of the type ^{1, exp^^), exp(2//!#)} etc. We see that minimizing Eq(|l0"D w.r.t. 
9 pushes up the energy as 



JU s a(0) 

enon + 2 i + (negligible terms) (11) 

\Hx\£ 



where 9* = l/(4|/ii|) log([/<j|/ii|£/4a(0)). The energy density continues to be inde- 
pendent of U s . Thus the growth of 9 stops on the scale of log(l/£), and we find 
a* = \\J a(0)U s C\iJLi\. Since a* does not have true LRO, our original assumption is 
validated, and hence we have a consistent solution. The solution at half filling thus has 
Quasi LRO since < A^A >~ £ 3//2 , and is superconducting in the sense of having power 
law order, as e.g. in a 2-d superconductor at finite temperatures. 

The d-wave superconductor: We now ask the question whether this class of wave 
functions works for d-wave superconducting correlations as well in 2-d. This case is clearly 
of topical interest in the context of the 2-d High T c cuprates. There is one possibility 
inherent in the discussion of Ref(|5|]), wherein we take B = J2T]kbk, with rjk = sign(wfc), 
and Wk = cos(k x ) — cos(k y ). This choice of the relative phases of the k state Cooper pairs 
preserves the important commutation relation [B, B>] = £ — N (where N is the number 
operator), since r/l = 1, and all the arguments of Ref(||) can be taken over, exactly as 
for the case of extended s-wave pairing discussed above. The conjugate variable that is 
enhanced now is A = —2J2£kVkbk, and the two particle wave function has the signs that 
are usually associated with d-wave pairing. The "gap" within a weak coupling approach 
in this case may be taken to be A = Aoe/^, and thus has jumps in k space where 
one expects nodes. Of course this may not be a fatal difficulty since there is no clear 
connection between the ground state obtained within the present framework, and the low 
lying excitation spectrum within a BCS framework. 

Within the variational framework, we can admit d-wave pairing provided the lattice 
has a nonzero orthorhombic distortion, i.e. the band dispersion is asymmetric in x and y 
directions. We can model this by setting = — 2(1 — e ) cos(k x ) — 2(l+e ) cos(k y ) with e a 
measure of the orthorhombicity. We choose a variational wave function exactly as in Eq(^), 
with S dw — S( w fci — w k 2 )b\bk 2 an d define the d-wave pairing operator D = zZw^bk- 
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With this we can essentially borrow all the results of the previous calculation provided 
we replace \i n — > is n , where v n = \ < T n > and T n = J^^^k- The fermi groundstate 
has the symmetry of the and hence we see the need for the distortion, for in its 
absence, v n = for all n > 1. We must also replace the initial conditions appropriately, 
thus A(#) = i < D^D > 9 , and A(0) = [wlfk]k- We need to examine the reality of 



A = y^i + (1 - p)v 2 - In case A is imaginary, the nature of solution changes drastically, 
and we have a limited enhancement, since the functions now oscillate as trigonometric 
functions, as 9 varies. If Ao is real, we can repeat everything said before and qualitatively 
find the same answers as for the extended s-wave case. The numerics require a knowledge 
of the elements v n and we will not discuss them further here. At half filling the state is 
similarly obtained, and the energy is different from that of the s-wave case only in terms 
of o{i/VZ), it is higher in general since v n are usually smaller than the \i n for n > 1. It 
must be remarked that in this case the correlation function a is also enhanced similarly 
to A, since these correlations get coupled. The correct symmetry of the order parameter 
is then of a mixed s and d type, rather than pure s or dQ, and indeed one should, in 
principle, optimize the form of the S operator by taking a suitable mixture of s and d 
functions to get the best possible energy. 

Summary and Discussion: In summary, we have presented variational wave func- 
tions that seem interesting on several counts. In many regards, we are tempted to say 
that the wave functions are at least as interesting as the models if not more. Firstly, 
the expectation values are computable analytically and throw light on an interesting and 
novel relationship, namely the enhancement in certain correlations at the expense of on 
site s-wave correlations. At half filling we have a state that is degenerate with the free 
fermi gas in the energy-density, and yet has critical superconducting correlations. Is this 
state the exact ground state of the model in Eq(]l|)? The coincidence of the energy density 
with the obvious lower bound e non , is not a clinching argument by itself, particularly in 
view of the main message of our work. The nature of the correlations, on the other hand, 
being of the form expected from the uncertainty principle arguments, provide a signifi- 
cant positive indicator. We feel that certain results found here are likely to be true, such 
as a* oc C 1 / 2 , and the form of energy corrections as in Eq ([TT|) , without however, being 
able to prove them in a rigorous fashion. Of the values of the coefficients of the leading 
behaviour, it is less easy to be certain. 

The critical behaviour of the Cooper correlation function < A^A > is of particular 
interest. It behaves as £2 at p = 1, and is oc otherwise, suggesting that the local 
density a{r) of the A operator has correlations < a^(r)a(0) >~ i^j^hir / 'cj(p)) , with a 
density dependent correlation length £(p) ~ ^ \ 2 /d , that diverges at half filling. Such a 
crossover is similar to that in the behaviour of the spin correlations in the Hubbard model 



in 1-d near half filling |T0[. Note that the behaviour of this function for the free fermi gas 
is ~ ^ , an d hence the correlations in this case are considerably shorter ranged than ours. 

We conclude by mentioning another nontrivial implication of the ideas discussed above, 
but going beyond the context of the models considered here. We have seen that the fermi 
gas on a lattice has the possibility of novel large enhancements in the ground state cor- 
relation functions of Cooper pairs, produced by repulsive terms of certain kinds. These 
are not reflected in either the total energy (and hence the compressibility) or the momen- 
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turn distribution function. The latter can be readily seen to be as sharp as in the ideal 
fermi case. Similarly one can check that the spin correlations are unaffected. The usual 
fermi liquid enhancements in spin and charge response that herald incipient instabilities 
are totally missing in these systems. As a result, any bosonic excitations available, in- 
cluding phonons, would find the electronic system much more prone to superconductivity 
of the usual sort than expected on grounds of the missing fermi liquid enhancements, 
and presumably be accompanied by enhanced T c 's. Estimates of such effects must await 
information on the spectral weight in the dynamical Cooper susceptibility. 
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Figure 1: The three correlations functions at p = 0.75 in 2-d vs. 6. 





Figure 3: The enhancement factor of the extended s-wave correlator a*/a(0) vs density 
in 2-d, and in 1-d (inset). 
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